Threshold model and the latest NA50 data on J/ip suppression in Pb+Pb collisions 
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Using the QGP motivated threshold model, where all the 
J/ip's are suppressed above a threshold density, we have an- 
alyzed the latest version of the NA50 data on the centrality 
dependence of the J/ip over Drell-Yan ratio. The data are 
not well explain in the model, unless the threshold density is 
largely smeared. Large smeared threshold density effectively 
excludes creation of any deconfined medium in the collision. 



I. INTRODUCTION 

In relativistic heavy ion collisions J / ip suppression has 
been recognized as an important tool to identify the pos- 
sible phase transition to quark-gluon plasma. Because 
of the large mass of the charm quarks, cc pairs are pro- 
duced on a short time scale. Their tight binding also 
makes them immune to final state interactions. Their 
evolution probes the state of matter in the early stage of 
the collisions. Matsui and Satz [1] predicted that in pres- 
ence of quark-gluon plasma (QGP), binding of cc pairs 
into a J/ip meson will be hindered, leading to the so 
called J/ip suppression in heavy ion collisions [1]. Over 
the years several groups have measured the J/ip yield in 
heavy ion collisions (for a review of the data and the in- 
terpretations see Refs. [2,3]). In brief, experimental data 
do show suppression. However this could be attributed 
to the conventional nuclear absorption, also present in 
pA collisions. 

In 1998 NA50 collaboration [4] published the results of 
centrality dependence of J/ip suppression in 158 GeV/c 
Pb+Pb collisions. Data gave the first indication of 
the anomalous mechanism of charmonium suppression, 
which goes beyond the conventional suppression in a nu- 
clear environment. The ratio of J /tp yield to that of Drell- 
Yan pairs decreases faster with Et in the most central 
collisions than in the less central ones. It has been sug- 
gested that the resulting pattern can be understood in 
a deconfinement scenario in terms of successive melting 
of charmonium bound states [4]. Later, the data were 
analyzed in a variety of models, with or without the as- 
sumption of deconfining phase transition [5-7,?, 8]. 

Recently, in Quark Matter 2002, NA50 collaboration 
presented preliminary analysis of the Et dependence 
of the J/ip over Drell-Yan ratio, obtained in the 2000 
Pb+Pb run [9]. The suppression pattern is changed, pre- 
sumably due to changed method of analysis. In 1998, 
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NA50 collaboration performed two types of analysis, the 
standard analysis where the J/ip and Drell-Yan cross sec- 
tions were measured directly to obtain the ratio. In the 
other analysis, called the minimum bias (MB) analysis, 
the Drell-Yan cross sections were replaced by the mini- 
mum bias cross section, to obtain, 
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The data beyond 100 GeV were obtained entirely from 
the MB analysis [4] . The latest NA50 data are obtained 
from the standard analyzed only. Compared to 1998 data 
[4], 2000 data [9] are flatter, suppression being more at 
low Et and less at high Et- However, the suppression 
obtained is still anomalous in the sense that normal nu- 
clear absorption model fails to explain it. In Quark Mat- 
ter 2002, NA50 collaboration also presented their anal- 
ysis of the nuclear absorption of J/ip in high statistics 
450 GeV pA collisions [10]. They estimated the J/ip nu- 
cleon absorption cross section (cr^lf N ) in the framework 
of Glauber model. High statistics 450 GeV pA data yield 
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— 4.4+1.0 mb [10]. They also estimate a common 

^ from latest pA and NA38 200 GeV/c S+U data 

[11], u^lf N =4.4 ± 0.5 mb. The extracted absorption 
cross section is much smaller than the earlier value of 6.4 
± 0.8 mb extracted from fit to earlier NA50 data [12] or 
7.1 ± 3.0 mb obtained from a fit to NA38 S+U data [11]. 
Within error, the S+U cross sections are compatible with 
pA cross sections. 

The changed suppression pattern in the latest NA50 
data [9], along with the high statistics pA data [10], with 
the implication of smaller J/f/j-nucleon absorption cross 
section, necessitated re-examination of the models, which 
were successful in explaining the earlier version of the 
NA50 data. Recently we have analyzed the latest NA50 
data in the QCD based nuclear absorption model [13]. 
The parameters of the model were fixed from the recent 
high statistics pA data. Without any free parameter the 
model could explain the latest NA50 data. Capella et al 
[14] also analyzed the data in the comover model. J/ip- 
nucleon absorption cross section was fixed at 4.5 mb. The 
latest NA50 data are explained with comover- J/ip ab- 
sorption cross section, <7 CO =0.65 mb. 

Blaizot et al [5] proposed the QGP motivated thresh- 
old model to explain the earlier version of the NA50 data 
[4] . To mimic the onset of deconfining phase transition 
above a critical energy density and subsequent melting 
of J/^'s, J/ip suppression was linked with the local en- 
ergy density. If the energy density at the point where 
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J/ip is formed, exceeds a critical value (e c ), J/ip's disap- 
pear. The critical energy density was then related with 
the (transverse) density and its value was obtained from 
a fit to the NA50 experimental data [4]. The NA50 data 
were well fitted with critical (threshold) density n c =3.7- 
3.75 fm~ 2 and <j^ N =6A mb. The model needs to be 
tested against the latest NA50 data [9], with changed 
suppression pattern, taking into consideration the new 
data on pA collisions [10]. 

Aim of the present paper is to analyze the latest NA50 
data [9] on the centrality dependence of J/ip over Drcll- 
Yan ratio, in the threshold model [5]. Plan of the paper 
is as follows: In section 2, we present a brief description 
of the threshold model. In section 3, NA50 data are 
analyzed in the model. Summary and conclusions are 
drawn in section 4. 



II. J/ip SUPPRESSION IN THRESHOLD MODEL 

The details of the threshold model could be found in 
[5]. In addition to the 'conventional' nuclear absorption, 
Blaizot et al. [5] introduced an anomalous suppression 
factor, S anom , such that all the J/ip's arc totally sup- 
pressed above a critical (threshold) density n c . The J/ip 
cross section at an impact parameter b as a function of 
Et is then written as, 
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S anom (h,s)P(b,E T ) 



where T e **(b) is the effective nuclear thickness, 

r-"(b) = j^dzp^z^pi-^ir r d ZIP {h,z>). 

For the density p, we use a Woods-Saxon form 



p(r) 



l + exp((r - R)/aY 



d 3 rp(r) = A (3) 



In [5], Blaizot et al used R = 1.1 A 1/3 = 6.52 fm and 
a=.53 fm. However, nuclear absorption has a sensitive 
(exponential) dependence on density and it is better to 
use more realistic values for the radius and the diffu- 
siveness parameters. For Pb we use, R = 6.624/to and 
a=0.549 fm [15]. c^lf N is the J/i/>-nucleon absorption 
cross section for which a value of =4.4 mb [10] is used. 

In Eq.2, P(b, Et) is the Et — b correlation function for 
which a Gaussian form is used [5], 

P(b, E T ) cx exp(-(E T - qN p (b)) 2 /2q 2 aN p (b)), (4) 

where N p (b) is the number of participant nucleons at 
impact parameter b. a and q are parameters related to 
dispersion and average transverse energy. For Pb+Pb 
collisions the parameters are, a— 1.27 and q=0.274 GeV 
respectively [5]. 

Sanom (b, s) in Eq.2 is the anomalous suppression fac- 
tor introduced by Blaizot et al [5]. They considered two 



forms for S anom . Assuming that all the J/^'s get sup- 
pressed above a threshold density (n c ), the anomalous 
suppression factor can be written as, 



>(b,s) = 9(n(b,s) - n c ) 



(5) 



where n is the transverse density (Eq.3 of [5]). In ref. 
[5] it was seen that if the theta function is smeared by 
a small amount, at the expense of another parameter, 
such that suppression is gradual rather than abrupt, the 
quality of fit to data improves considerably. This was 
done by writing, 



( (b, s) = 0.5[1 - tanh(X(n(b, s) - n c ))] (6) 



In both the forms, effect of Et fluctuations at a fixed 
impact parameter is taken into account by rescaling the 
density as, n — > tiEt/ < Et > (b) [5]. Threshold model 
parameters n c and A can be obtained by fitting the latest 
NA50 data on centrality dependence of J/tjj over Drell- 
Yan ratio. 

The Drell-Yan pairs do not suffer any final state inter- 
action and the cross section at impact parameter b as a 
function of Et can be calculated as, 
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FIG. 1. The NA38 data on the centrality dependence 
of J /ip over Drell-Yan ratio in 200 GeV/c S+U collisions. 
The line is the Glauber model of nuclear absorption with 
a aba' N '=4.4 mb. For the Et — B correlation, we have used 
Gaussian form Eq.4, with a=3.2 and q=0.74 GeV [2]. 



III. COMPARISON WITH LATEST NA50 DATA 

For comparison with NA50 experimental data on the 
centrality dependence of J/ip over Drell-Yan ratio, nor- 
malizing factor N = B^^a'^ /a^Jj is needed. It can 
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also be considered as a fitting parameter and obtained 
along with the threshold model parameters, n c and A. 
However, normalizing factor and threshold model param- 
eters are correlated. Higher normalizing factor can be 
compensated by increasing the threshold density. Not- 
ing that the ratio of Pb+Pb to S+U normalizations is 
equal to 1.051 ± 0.026 [10,14], NA38 data on the cen- 
trality dependence of J/tp over Drell-Yan ratio in 200 
GeV S+U data [11] can be fitted to obtain the normal- 
izing factor for S+U collisions and can be rescaled to 
obtain the same for Pb+Pb collisions. In S+U collision, 
dcconfincment transition is not expected. The observed 
suppression is due to nuclear absorption only. Fixing 

a^[f N =4A mb, NA38 data are fitted in the conventional 
nuclear absorption model to obtain the normalizing fac- 
tor. In Fig.l, NA38 data along with the best fit obtained 
with B^a^x /(Jhh=39.65 is shown. Within the error 
bars, the data are well explained in the nuclear absorp- 
tion model. Normalizing factor for Pb+Pb collision is 
then obtained by scaling the normalizing factor for S+U 
collisions by 1.051 [10]. 

We may mention that recently, Capella et al [14] an- 
alyzed the S+U data in the comover model. With 
a abt^ '=4-5 mb, data allow for a small comover interac- 
tion, <7 CO =0.65. For the B^^a^/^ /a^Jj they used a value 
of 47, obtained from a fit to latest NA50 Pb+Pb data 
in the comover model. However, as we find the data are 
well explained in terms of nuclear absorption only, use of 
comover interaction in S+U collision is debatable. 
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FIG. 2. The latest NA50 data on the centrality depen- 
dence of J ftp over Drell-Yan ratio in Pb+Pb collisions. The 
dotted line is the Glauber model of nuclear absorption with 
a J l^ N '=4.4 mb. The dash-dotted line is the ratio in the 

abs 

threshold model with n c =3.61 fm~ 2 . The dash-dot-dot line 
is obtained with n c =3.72 fmT 2 and \—2 fm 2 (fixed). The 
solid line is obtained with , n c =3.82 fmT 2 and A=0.77 fm 2 . 

In Fig. 2, latest NA50 data [9] on the centrality depen- 



dence of J/tp over Drell-Yan ratio in Pb+Pb collisions, is 
shown. Just to show that the latest data are also anoma- 
lous, we have shown the Glauber model calculation with 
a abs'^ '=4-4 mb (the dashed line). Only for very periph- 
eral collisions, the Glauber model of nuclear absorption 
fits the data. For more central collisions, it produces 
much less suppression than the data exhibit. In Fig.l, 
the dash-dotted line is the best fit obtained to the data in 
the threshold model, without any smearing of the thresh- 
old density. We obtained threshold density n c =3.61 ± 
0.06/m~ 2 . The value is less than value n c =3.7 /m~ 2 , 
obtained by Blaizot et al [5] from the earlier version of 
the data, mainly due to lesser value of the J/V'-Nucleon 
absorption cross section. To produce similar suppression, 
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being less, anomalous suppression has to increase. 



Interestingly, we find that the threshold model with a 
single parameter does not give a proper description to 
the data. In the intermediate range of Et, agreement 
with data is not good. Smearing the threshold density 
by a small amount do not improve the quality of fit. In 
Fig. 2, the dash-dot-dot line is obtained with smearing 
the threshold density by a small amount (A fixed at 2 
fm 2 ). The best fit to the data is then obtained with 
threshold density,n c =3.72 fm~ 2 . Again the data are not 
well explained. In Fig. 2, the solid line is the best fit 
obtained to the data varying both the parameters. We 
obtain, n c =3.82 ± 0.09 fm,- 2 and A=0.77 ± 0.14 fm 2 . 
The model then reproduces the data through out the 
Et range. Small value of A required for good fit to data 
indicate that considerable smearing of the threshold den- 
sity is required for proper description of the NA50 data. 
The anomalous suppression is not abrupt but increases 
gradually with density. The analysis suggest that with 
a abs' N '=4.4 mb, unless the threshold density is largely 
smeared, threshold model do not give a proper descrip- 
tion of the latest NA50 data on the centrality dependence 
of the J/tp suppression. 



IV. SUMMARY AND CONCLUSIONS 

To summarize, the latest NA50 data on the central- 
ity dependence of J/tp over Drell-Yan ratio in Pb+Pb 
collisions are analyzed in the QGP motivated threshold 
model. In the threshold model, in addition to the con- 
ventional nuclear absorption, an anomalous suppression 
is introduced, such that above a threshold density, all 
the J ftp's are absorbed. To be consistent with latest pA 
data on J/tp absorption, we have used c^f^ '=4-4 mb. 
Threshold model with a single parameter, the threshold 
density (n c ) donot give a proper description of the cen- 
trality dependence of J ftp over Drell-Yan ratio. The best 
fit to data is obtained with threshold density n c =3.61 
fm~ 2 fails to explain the data in the intermediate Et 
range. If the threshold density is smeared at the expense 
of another parameter (A) the model could explain the 
latest NA50 data with n c =3.82 fm~ 2 and A=0.77 fm 2 . 
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Small value of X—0.77 fm 2 indicate that for a proper 
description to the NA50 data, the threshold density has 
to be smeared considerably. Thus onset of anomalous 
suppression is not sudden, rather gradual. Over a den- 
sity range of 2.4-5.2 /m~ 2 , anomalous suppression factor 
change from 0.9 to 0.1. Originally, threshold model was 
devised to mimic the melting of J ftp's in a deconfming 
medium. With nominal smearing of the threshold den- 
sity, the essence of the model is not lost. Large smear- 
ing of the threshold density, as required to fit the latest 
NA50 data, effectively excludes formation of deconfming 
medium. The medium where the anomalous suppression 
is taking place, do not melt the J/ip suddenly, but rather 
gradually, more like in a nuclear/comover environment. 
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